Interplay of density and phase fluctuations in ultracold one-dimensional Bose gases 
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The relative importance of density and phase fluctuations in ultracold one-dimensional atomic 
Bose gases is investigated. By defining appropriate characteristic temperatures for their respective 
onset, a broad experimental regime is found, where density fluctuations set in at a lower temperature 
than phase fluctuations. This is in stark contrast to the usual experimental regime explored up to 
now, in which phase fluctuations are largely decoupled from density fluctuations, a regime also 
recovered in this work as a limiting case. Observation of the novel regime of dominant density 
fluctuations is shown to be well within current experimental capabilities for both 23 Na and 87 Kb, 
requiring relatively low temperatures, small atom numbers and moderate aspect ratios. 
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Coherence plays a crucial role in the behavior of a sys- 
tem, and is central to our understanding of laser, matter 
wave and condensed-matter physics. The achievement of 
quantum degeneracy in ultracold one-dimensional (ID) 
atomic clouds either in magnetic 0, or optical traps 0, 
or maintained close to microfabricated surfaces known 
as 'atom chips' 3] has led to a unique experimentally 
controlable system for fundamental studies of coher- 
ence in low dimensional systems. Potential applications 
of such systems include 'on-chip' atom interferometers 
0, 0|, and continuously-operating atom lasers Q- ID 
systems differ from their 3D counterparts, in that low- 
dimensional systems are generally prone to large phase 
fluctuations, which tend to destroy long-range phase co- 
herence 0, H, 0, 0, El 0| • Nonetheless, harmonically- 
confined ID gases can maintain significant coherence 
across the system size at sufficiently low temperatures, 
whereas, at higher temperatures, coherence is limited to 
smaller regimes. In this case, the system contains a quasi- 
condensate Q. 

The suppression of coherence with reduced dimension- 
ality has been observed experimentally in very elongated 
3D Bose-Einstein Condensates (BECs) [3 El El [3 
17]. In these systems, the equilibrium coherence proper- 
ties were found to be largely insensitive to density fluctu- 
ations, even at temperatures close to the critical temper- 
ature, where most of the atoms are in the thermal cloud. 
This can only be true, provided phase fluctuations set 
in at temperatures much lower than density fluctuations, 
which is typically considered as the only attainable ex- 
perimental limit in such systems. 

In this paper, we demonstrate that the opposite 
regime, in which density fluctuations set in at a lower 
temperature than phase fluctuations, is actually within 
current experimental reach. We determine optimum con- 
ditions for reaching this novel regime for two atomic 
species, highlighting in particular the case of 23 Na, which 
appears to be a better candidate for such experiments. 
The interplay between density and phase fluctuations is 
examined by comparing appropriately defined character- 



istic temperatures for their respective onset. Such an 
approach has already been used in other systems to ex- 
amine the relative interplay of phase and density fluctu- 
ations of the order parameter of the system. A striking 
such example is the case of high-T c superconductors, for 
which the interplay between amplitude and phase fluctu- 
ations is believed to be related to the 'pseudogap' phase 
|l8j| . The present study probes the entire 'crossover' be- 
tween the previously observed regime of dominant phase 
fluctuations, and the currently unexplored regime with 
dominant density fluctuations in weakly-interacting ID 
ultracold atomic Bose gases. This is performed by vary- 
ing the number of atoms in the system (Fig. 1), in close 
analogy to the well-known phase-amplitude-fluctuations 
interplay studied in high-T c superconductors as a func- 
tion of doping EH • 
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FIG. 1: Ratio T^/T„ highlighting interplay between density 
and phase fluctuations, as a function of atom number. Line 
Ttf, — T n indicates crossover between regimes of dominant 
phase (Tcj, < T„) and density fluctuations. Circles: 23 Na re- 
sults for lu z = 2tt x 10 Hz, lu± = 2tt x 2500 Hz. 87 Rb results 
also shown for same trap configuration (dashed line) and for 
uj z — 2-7rx50Hz (squares) Inset: Corresponding characteristic 
temperatures (open symbols) and T n (filled). 
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In atomic gases, such characteristic temperatures can 
be determined as follows: The temperature for phase 
fluctuations, is obtained by measuring the normal- 
ized off-diagonal one-body density matrix, also known in 
the language of quantum optics as the first order cor- 
relation function. Roughly speaking, corresponds to 
the temperature at which the latter quantity decays to 
zero at the edge of the quasi-condensate Q , thus marking 
an approximate boundary between the regimes of 'true' 
condensation and quasi-condensation [l(l |. On the other 
hand, the temperature for density fluctuations can be 
determined by measuring the quasi-condensate fraction. 
In particular, we define T„ as the temperature at which 
there is a 10% quasi-condensate depletion, since this can 
be readily measured experimentally. For a given experi- 
mental configuration, both and T n , and therefore their 
ratio, are fixed, and experiments so far have only been 
performed in the limiting regime T^/T n -C 1 (see later). 
However, the flexibility of atomic gas experiments en- 
ables complete tunability of the ratio T < j > /T ni thus al- 
lowing one to probe the full interplay between density 
and phase fluctuations. This is explicitly shown in Fig. 
1, for the particular example of variable atom number 
at fixed trap configuration. As evident from this figure, 
the novel regime of dominant density fluctuations, which 
can be approximately defined by > T n , is attainable 
for a broad range of experimental parameters, both for 
23 Na and 87 Rb. The corresponding characteristic tem- 
peratures for both species are shown in the inset. 

The equilibrium properties of the system are expected 
to depend on the interplay between and T n . In par- 
ticular, the role of density fluctuations on the coherence 
of the gas can be determined by looking at the difference 
in the measured coherence length of a depleted quasi- 
condensate, from calculations based on the correspond- 
ing undepleted system with the same total atom num- 
ber. Experiments to date have inferred a negligible ef- 
fect of density fluctuations on the equilibrium coherence 
properties 0, Q, 0, 0] , based on good agreement be- 
tween measured data and a theory which neglects den- 
sity fluctuations 0, |2l|. Surprisingly, this was argued 
to be the case even at temperatures close to the effec- 
tive 'transition' temperature, where most of the quasi- 
condensate is known to be depleted by thermal excita- 
tions. In this paper, we demonstrate that such a con- 
clusion is only correct in the regime -C T n , which is 
in fact the only regime explored experimentally to date. 
Even within this regime, more recent work [l7| indicates a 
currently unaccounted for reduction of 20% in the coher- 
ence length, which we attribute to density fluctuations. 
We also demonstrate that this reduction can be largely 
enhanced by appropriate choice of parameters. 

Theory: Early theoretical work addressing trapped ID 
systems treated phase fluctuations exactly, but neglected 
density fluctuations 0, thus applicable only in the limit 
of an undepleted quasi-condensate. Quasi-condensate 
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FIG. 2: Top: Density profiles at T — T$ for the opposite 
regimes of T^/Tn — (a) 0.5 and (b) 2.1: Quasi-condensate 
(solid black lines) and thermal (dashed black) density profiles 
and corresponding zero-temperature profiles (grey) at fixed 
total atom number (N = (a) 3500, and (b) 300). Densities 
scaled to the zero-temperature central density and distance 
scaled to the zero-temperature Thomas- Fermi radius Rtf(0). 
Bottom: Corresponding normalized first-order spatial corre- 
lation functions g^'(0, z/Rtf(0)), evaluated at the trap cen- 
tre, including (black) and excluding (grey) quasi-condensate 
depletion. Plotted data correspond to the highlighted 23 Na 
points of Fig. 1, with R T f{T)/R T f(0) = (a,c) 0.9 and (b,d) 
0.7 and quasi-condensate fractions Nq/N — (a) 0.96, (b) 0.73. 



depletion was subsequently included in a self-consistent 
treatment of both phase and density fluctuations, lead- 
ing us to the construction of the ID phase diagram in 
the weakly-interacting limit |10j . The resulting theory 
was shown to be free of both infrared and ultraviolet di- 
vergences and therefore valid in any dimension. Impor- 
tantly, and in stark contrast to other existing theories 
Ell ll^ - such an equation of state enables a direct 
ab initio determination of quasi-condensate density pro- 
files. In the harmonic trap, densities are obtained in the 
local density approximation, the effective temperature- 
dependent quasi-condensate size Rtf(T) is obtained in 
the Thomas-Fermi limit, and the coherence of the system 
is investigated via the first-order normalized correlation 
function. 

A ll q uasi-condensate experiments performed to date 
[l3L Il4 Hil IT 5 Il7i have been analyzed on the basis of 
the theory of [!J,|21| which entirely neglects both quantum 
and thermal fluctuations. Although the effect of quan- 
tum fluctuations is small for experimentally relevant con- 
ditions, thermal fluctuations clearly become important 
with increasing temperature, as evident from the exper- 
imentally observable quasi-condensate depletion. It is 
therefore not a priori clear why such density fluctuations 
have not been observed to couple to the phase fluctua- 
tions of the system. In order to ascertain the importance 
of density fluctuations on the equilibrium properties of 
the system in an ab initio manner, we henceforth compare 



T„ / T, 



the theory of upon neglecting the quantum fluctu- 
ations, to the conventional ID theory which additionally 
ignores thermal quasi-condensate depletion 0, 0] . The 
general expression for the first-order correlation function 



exp 



(-([x(z)-x(0)f>/2 



thus becomes g^'(0,z) 

where x{z) i s the phase operator, satisfying |lf 



(KW-Je(o)n = 
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Here Pj{z) are Legendre polynomials of order j, and 
A j = v/fjIl/V/b, Bj = (y/(j + l/2)hLo 3 /J7)/2. 
The frequencies are given by Uj = + l)/2 w z , /i' 

is the 'renormalized' chemical potential including quasi- 
condensate depletion, k the effective ID scattering length 
and l z the harmonic oscillator length corresponding to 
longitudinal confinement lj z . N(Tiojj) is the usual Bose 
distribution function. Density fluctuations are explicitely 
maintained in the above expression, via the temperature- 
dependent quasi-condensate size Rtf{T) appearing in 
the prefactor. The general expression quoted above can 
be readily reduced to the conventional theory 0, [nj 
which ignores quasi-condensate depletion by replacing 
Rtf(T) by the corresponding zero temperature quasi- 
condensate size -Rtf(O) at the same total atom number, 
and ignoring the B 3 contributions. Since Rtf(T) oc 
the former step is equivalent to replacing the renor- 
malized chemical potential at temperature T by the 
corresponding zero-temperature one for the same to- 
tal atom number. Then, the 'classical' approximation 
N{TiLOj) w ksT /hujj leads to the definition of the char- 
acteristic temperature = (hui z ) 2 N/kBfJ- 0- All pre- 
sented results maintain the full Bose distribution func- 
tion N(hujj) — [exp(/3hujj) — 1] , and are in the weakly- 
interacting regime, with 7 = l/(n£) 2 < 10~ 3 

Results: The effect of density fluctuations on the equi- 
librium properties of the system is best assessed by in- 
vestigating the latter at the characteristic temperature, 
T$, where phase fluctuations set in. This study, shown 
in Fig. 2, reveals marked differences in both density pro- 
files (top figures) and correlation functions (bottom) be- 
tween the regimes of dominant phase and density fluc- 
tuations. In particular, inclusion of density fluctuations 
leads to significant quasi-condensate depletion, reducing 
both the central quasi-condensate density and the quasi- 
condensate size Rtf(T) < Rtf(0)- Moreover, within 
the region of validity of the finite temperature correla- 
tion functions, i.e., z < Rtf(T), density fluctuations are 
shown to lead to a potentially significant reduction in the 
coherence. 




FIG. 3: (a) Fractional deviation &L coh (T) = (L° coh (T) - 
LZ h (T))/L? oh (T) of the coherence length L° coh (T) ignoring 
density fluctuations, from actual result L^ oh (T) including den- 
sity fluctuations, plotted as a function of the 'crossover' pa- 
rameter T^/Tn (bottom axis) and T^/T c (top). Curves eval- 
uated at T — Tij, (open triangles) and T — 2T$ (closed), (b) 
Dependence of 5L co h{T) on scaled temperature T/T n (bot- 
tom), or equivalently T/T c (top). Parameters as in Fig. 2. 



To quantify this reduction further, Fig. 3 shows the 
dependence of the fractional deviation of the coherence 
length induced by the presence of density fluctuations 
(a) on the 'crossover' parameter T^/Tn, and (b) on ap- 
propriately scaled temperature. The coherence length is 
defined as the distance at which the correlation function 
decays to half its value, i.e. (0, L co h(T)) = 0.5. The 
neglect of density fluctuations leads to a large overesti- 
mate of the coherence length. Importantly, this occurs 
already at the crossover region « T n even at tem- 
peratures as low as at which density fluctuations are 
usually considered negligible Using the approximate 
scaling T n w 0.2T C valid for the examined parameters, we 
conclude from Fig. 3(a) that density fluctuations should 
be observable at very low temperatures, even within ex- 
perimental uncertainties, provided T^/T c > 0.1. The ap- 
proximately linear increase in temperature of 5L co h(T) 
plotted in Fig. 3(b) is universal when temperatures are 
scaled in terms of T n , or equivalently, T c . However, cor- 
responding features differ noticeably in terms of T/T^,, 
suggesting that is not the most appropriate scaling 
parameter for such studies. 

The dependence of the coherence length, scaled to the 
zero-temperature system size, on appropriately scaled 
temperature is shown in Fig. 4. Filled (open) symbols 
correspond to the cases T$ > T n (T^ < T n ), with (black) 
or without (grey) density fluctuations included. The 
coherence length decreases approximately exponentially 
with temperature, consistent with experimental findings 
[lfij . Importantly, however, in both cases, density fluctu- 
ations lead to a considerable shift of the curves to lower 
values. This is also true when temperatures are scaled 
in terms of (Fig. 4(b)), as usually done experimen- 
tally 0, Q, 0, Q, E3 j m which case results without 
density fluctuations fall roughly on a universal curve (in- 
set). However, inclusion of density fluctuations is found 
to lead to a noticeable downward shift of the correlation 
function, already in the regime < T n /2 (black cir- 
cles). The magnitude of this shift clearly increases with 
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FIG. 4: Temperature dependence of the scaled coherence 
length L coh {T)/R T F(0) on (a) T/T c and (b) T/T for param- 
eters of Fig. 2: T^/Tc « 0.1 (open symbols) and 0.4 (filled), 
without (grey) or with (black) density fluctuations. Inset: 
Deviation in regime 0.4 < T/T$ < 1.4. 

increasing T^/T„, as demonstrated by the black squares. 
Note that results with density fluctuations are fully con- 
sistent with corresponding results based on a stochastic 
Langevin approach in the classical field approximation 
|2(i|. to be presented elsewhere. 

Current Experiments: Measurements of phase fluctu- 
ations to date have only been performed close to the 
ID regime. In these experiments, the departure from 
pure ID excitations arising from /i ~ few Tko±_ and 
the small values of the phase degeneracy temperature 
(T^/Tc <C 0.1) tend to suppress density fluctuations. 
Early experiments based on interferometric techniques 
(lif . or Bragg spectroscopy ^(| found no evidence of 
density fluctuations on the equilibrium coherence proper- 
ties, within their observational uncertainties. Nonethe- 
less, more recent measurements of the temperature de- 
pendence of the scaled coherence length performed at 
temperatures T/T^ < 6 17] cannot be interpreted by a 
3D theory which ignores density fluctuations [2l| unless 
the corresponding predictions are shifted to lower val- 
ues by about 20%. We believe that this shift arises, at 
least partly, from the effect of density fluctuations, as 
supported qualitatively by our present ID analysis. 

Proposed Experiment: The crossover between the 
regimes of dominant phase and density fluctuations can 
be probed in a controlled manner by repeated measure- 
ments of the equilibrium coherence properties of the sys- 
tem at fixed temperature and variable atom number. In- 
stead of decreasing the atom number, one could achieve 
the same effect by reducing the scattering length via 
a Feshbach resonance |22j, since both these effects de- 
crease the system nonlinearity. The 'crossover' param- 
eter determining the relative interplay between phase 
and density fluctuations was numerically found to obey 
T<p/T n w r](hLu z /fi). We first focus our discussion on 
23 Na, which appears to be a good candidate for such 
experiments, and for which r\ sa 25. Based on the 
usual 3D scattering length of 2.75nm, this gives T^jT n w 
10 4 Hz 2 / 3 x (w*/wi)V3/JV 2 /3. 



The proposed regime is achieved when a number of 
competing conditions are satisfied: Firstly, ID thermal 
excitations require kT < huj±, implying a relatively large 
transverse confinement of the order of a few kHz, which is 
easily achievable in experiments with atom chips Q . The 
condition for ID quasi-condensation, (i < fko±, translates 
into a condition for the maximum number of atoms at 
given trap confinement, namely N 2 : (w 2 / 'u j_) < C\. Fi- 
nally, the requirement for dominant density fluctuations, 
T$ > T n , places a competing condition N (iv^/uJz) < C%. 
For 23 Na, we find C\ ~ few x 10 8 Hz and C2 ~ few x 
10 12 Hz. At the chosen confinement lu z = 2tt x 10Hz and 
uj± = 2ir x 2.5kHz, the ID quasi-condensate condition is 
the most restrictive for the maximum number of atoms 
required, whereas the opposite is true for the same aspect 
ratio uj±/lo z but larger oj±_ = 2ir x 10kHz. 

Throughout this work, we have focused on 23 Na, on 
the grounds that such experiments may be easier because 
23 Na enters the regime of dominant density fluctuations 
with a relatively large atom number N > 1000 for a 
broad parameter range. However, this comes at the ex- 
pense of a stringent requirement on low temperatures, 
typically around 20nK. The same regime is also attain- 
able with 87 Rb, for which 77 » 10. Although 87 Rb re- 
quires lower aspect ratios and higher temperatures than 
23 Na, it nonetheless requires a reduced number of atoms, 
with N < 700 even for optimized aspect ratios. The op- 
timum candidate for such an experiment thus depends 
on the flexibility of the experimental set-up, with respect 
to cooling efficiency and atom detection techniques. 

The regime of dominant density fluctuations typi- 
cally requires the phase fluctuation temperature, T^, to 
be comparable to the effective 'transition' temperature 
T c . In this regime, the suppressed phase fluctuations 
are dominated by their coupling to density fluctuations. 
Such a regime may prove useful for potential applications, 
due to the presence of large coherence even at relatively 
high temperatures T ~ 0.5T C . However, the effect of 
density fluctuations on the coherence length can already 
be observed at the much weaker approximate condition 
Tff, > 0.1T„, which increases the value C2 by a factor of 
10. 

In conclusion, we proposed a new regime of quasi-one- 
dimensional weakly-interacting atomic gases, in which 
density fluctuations set in at a lower temperature than 
phase fluctuations. This regime was shown to be experi- 
mentally accessible for small atom numbers, low temper- 
atures and moderate aspect ratios. It can be experimen- 
tally identified by the noticeable reduction in the equi- 
librium coherence length of the quasi-condensate, whose 
extent depends critically on the ratio of the phase degen- 
eracy temperature to the effective transition tempera- 
ture T c . Given the pronounced nature of density fluctua- 
tions in homogeneous systems, an alternative observation 
of dominant density fluctuations can be performed in the 
new box- like trap |23| . The unprecedented experimental 
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control offered by ultracold atomic gases may offer fur- 
ther insight into the interplay between phase and density 
fluctuations in analogous condensed-matter systems. 
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